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^ ! Abstract 

\o : 

C\l . We consider DM that only couples to SM gauge bosons and fills one gauge multiplet, 

e.g. a fermion 5-plet (which is automatically stable), or a wino-like 3-plet. We revisit 
r£3 • the computation of the cosmological relic abundance including non-perturbative 

corrections. The predicted mass of e.g. the 5-plet increases from 4.4 TeV to 10 TeV, 
Oh| and indirect detection rates are enhanced by 2 orders of magnitude. Next, we 

c-| \ show that due to the quasi- degeneracy among neutral and charged components of 

the DM multiplet, a significant fraction of DM with energy E> 10 17 eV (possibly 



O ; 

^ ; Contents 

o : 

t^- . 1 Introduction 2 

o ■ 

>• ■ 2 Non perturbative Sommerfeld corrections 2 

^ ! 2.1 Computing Sommerfeld corrections 5 

& : 

3 Cosmological abundances 7 

4 Indirect detection of Minimal Dark Matter 13 

5 Minimal Dark Matter at Ultra High Energies 15 

5.1 Propagation of UHE Minimal Dark Matter 15 

5.2 Production of UHE Minimal Dark Matter 18 

5.3 Detection of UHE DM ± 20 

6 Conclusions 20 
A Annihilation cross sections 22 



1 



1 Introduction 



Observations and experiments tell that: 

i) Dark Matter (DM) exists [1], with abundance n DM h 2 = 0.110 ± 0.005 [2]; 

ii) DM is cosmologically stable; it has no electric charge [3], no strong interactions [4] and 
the DM coupling to the Z boson is smaller than ~ 10~ 3 (?2 [5]. 

i) suggests that DM might be a weak-scale particle: indeed fluM ~ 1 is the typical freeze-out 
abundance of a typical particle with weak-scale mass M and coupling g: M/g ~ ^T Q Mp\ ~ 
TeV. ii) can be satisfied in appropriate models, although these properties are not typical weak 
scale particles. 

We here want to study Minimal DM (MDM) models: we add to the Standard Model 
(SM) one new multiplet, that only has gauge interactions and a SU(2)L-invariant mass term 
M. In this context, ii) singles out scalars or fermions with zero hypercharge that fill a SU(2)l 
representation with odd dimension: n = 3, 5, ... [6]. The neutral DM component is then lighter 
than the charged DM ± components, by only 166 MeV. DM interacts with ordinary matter 
via one-loop exchange of W ± : future direct DM searches will hopefully reach the sensitivity 
needed to detect it. 

The fermion 5-plet is particularly interesting: all interactions other than the gauge inter- 
actions are incompatible with gauge and Lorentz symmetries [6], so that it is automatically 
stable thanks to an accidental symmetry, like the proton in the SM. Other candidates that we 
will consider lack this feature (i.e. some extra symmetry is responsible for their stability) but 
are also interesting. The fermion 3-plet is well known from supersymmetry with matter-parity, 
as 'pure wino' DM; the fermion doublet case as 'pure higgsino' DM. We will also study scalar 
MDM candidates. In general, this scenario has one free parameter, the DM mass M, already 
fixed by fi DM , such that all DM signals are univocally predicted. We here address the following 
issues. 

First, we complete the computation of the cosmological abundance performed in [6] taking 
into account p-wave annihilations and renormalization of the gauge couplings up to the DM 
scale, and, more importantly, we include non-perturbative Sommerfeld corrections: the DM 
wave-functions get distorted by Coulomb-like forces mediated by SM vectors. Their relevance 
was pointed out in [7], in the case of gluino and wino as lightest supersymmetric particle. 
In section 2 we discuss the basic physics of non-perturbative Sommerfeld corrections, and in 
section 3 we show our results for the cosmological DM abundance. 

Next, in section 4 we compute the DM°DM° annihilation rates relevant for indirect MDM 
signals, including the sizable enhancement due to Sommerfeld corrections. 

Finally, in section 5 we study the possibility that the ultra-high-energy cosmic rays (UHE 
CR) contain some DM particle (although this looks difficult within the standard CR acceleration 
mechanism), showing that MDM candidates can cross the earth arriving to the detector in the 
charged state. We discuss how it could manifest in detectors such as neutrino telescopes. 

Section 6 presents our conclusions. 

2 Non perturbative Sommerfeld corrections 

If the DM mass is M m Mz and if the DM coupling is g ~ g2, DM DM annihilations into 
SM vectors have a too large cross section, and consequently give a too low freeze-out DM 
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abundance. Since f^DM 8S M for M > Mz, the observed DM abundance is obtained for a value 
of M sufficiently larger than Mz- 1 Since this value turns out to be in the TeV or multi-TeV 
range, we can ignore SM particle masses when computing the annihilation rates. Using SU(2) 
algebra, ref. [6] obtained a single expression for the cosmological abundance of any MDM 
candidate, taking into account all (co-)annihilations in s-wave approximation. We here add 
p-wave annihilations and renormalization of the gauge couplings up to the DM mass M: each 
one of these effects gives a (9(5%) correction to Qbm- Precise formulae are given in appendix A. 

Furthermore, we take into account non-perturbative electroweak Sommerfeld corrections. 
Their relevance was pointed out in [7] and might look surprising, so we start with a general 
semi-quantitative discussion before proceeding with the detailed computations. Scatterings 
among charged particles due to point-like interactions are distorted by the Coulomb force, 
when the kinetic energy is low enough that the electrostatic potential energy is relevant. This 
leads e.g. to significant enhancements of the /i~/U. + annihilation cross section (attractive force) 
or to significant suppressions of various nuclear processes (repulsive force). These effects can be 
computed with a formalism developed by Sommerfeld [9]. In the language of Feynman graphs, 
these effects are described by multi-loop photon ladder diagrams. Perturbative computations 
would include only the first few diagrams: the Sommerfeld enhancement is non-perturbative 
because a resummation of all ladder diagrams is needed. 2 The generalization of the Sommerfeld 
formalism to the case of DM DM annihilations, that involves non abelian massive vectors, was 
presented in [10]. 

Let us first discuss Sommerfeld corrections due to one abelian vector with mass My and 
gauge coupling a: this case already contains the relevant physics and can be analyzed in a 
simpler way. Cross sections at low energies are dominated by s-wave scattering, with p-wave 
giving corrections of relative order T/M. Since the DM DM annihilation is local (i.e. it occurs 
when the distance between the two DM particles is r ~ 0), the non perturbative correction is 
given by R = |?/>(oo)/?/>(0)| 2 , where ip{r) is the (reduced) s-wave-function for the two-body DM 
DM state with energy K, that in the non-relativistic limit satisfies the Schrodinger equation 

1 ^ + V ■ ij> = V = ±-e~ M - r (1) 



M dr 2 r 

with outgoing boundary condition ip f (oo)/ip(oo) ~ iM/3. Here K = M/3 2 is the kinetic energy 
of the two DM particles in the center-of-mass frame, where each DM particle has velocity (3. 
The — sign corresponds to an attractive potential and the + sign to a repulsive potential. We 
define e = My/M, the adimensional ratio between the vector mass and the DM mass. 

Eq. (1) is the prototype of the equations that describe corrections to DM°DM° annihilations 
mediated by a Z boson, or DM + DM _ or DM + DM + co-annihilations mediated by a 7 or W 
or Z, or corrections mediated by a gluon (in models where coannihilations with and between 
colored particles are relevant). In the three cases, the coupling a that appears in eq. (1) would 
be «2, <y C m and a s respectively. In the non abelian case, eq. (1) becomes a matrix equation and 



1 Another solution exists for M sufficiently lower than Mz, such that annihilation into vectors are kincmati- 
cally suppressed. In the minimal scenario that we consider, this solution is excluded by LEP2 and other collider 
data, but in general it is still allowed in some part of the parameter range of non- minimal scenarios [8] . 

A DM mass M « Mz can of course also be obtained by reducing g, i.e. by assuming that DM is any SU(2)l 
multiplct appropriately mixed with a singlet. 

2 Sometimes in the literature 'non-perturbative' is used with a different meaning: to denote effects that 
vanish in the perturbative limit because suppressed e.g. by e _1 / Q factors. 
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Figure 1: Iso- contours of the non-perturbative Sommerfeld correction to the DM DM anni- 
hilation. Here a is the coupling constant, (3 is the DM velocity, e is the ratio between the 
vector mass and the DM mass. The labels indicate where some classes of DM candidates lie 
in this plane: 'weak 7 indicates weak-scale DM particles, 'TeV indicates DM with multi-TeV 
mass, and 'strong' indicates strongly-interacting particles that in some models give dominant 
co- annihilations . Within the shaded region thermal masses dominate over masses, effectively 
shifting the value of ct/e as indicated by the arrow. 



V is the sum of the various contributions mediated by all SM vectors. Higgs exchange can also 
be relevant, in models where DM sizably couples to the Higgs. 

For e = (massless vector) the Schrodinger equation has the same form as the one that 
describes e.g. the hydrogen atom, and it can be analytically solved: the Sommerfeld factor R 
that multiplies the perturbative cross section is 

R = Y —^ X x = ±- (forM y = 0) (2) 

This shows that R sizably differs from 1 at (3 < not. DM DM annihilations into SM particles 
freeze-out when the temperature T cools down below T/M ~ l/ln(Mpi/M) ~ 1/26. This 
happens to be numerically comparable to the SM gauge couplings, i.e. c*2 ~ 1/30. Consequently 
DM freeze-out occurs when (3 ~ 0.2 i.e. ttx ~ 1: the Sommerfeld correction is significant, 
R-O(l). 

Of course we must take into account that the relevant W, Z vectors are massive: since R 
must now be computed numerically, it is convenient to first identify on which parameters R 
depends. We notice that R only depends on the two ratios a/ (3 and a/e, so that R can be 
plotted on a plane. This can be proofed noticing that R is adimensional and that physics is 
invariant under r — > Ar, M M/X, My — > Myj\. Fig. 1 shows iso-contours of R as function 
of at/fl and a/e. We can distinguish various regions: 
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a) Non-perturbative effects are negligible (i.e. \R — 1| -C 1) if 'kol <C and/or if a <C e. We 
have seen that na ~ f3, so non-perturbative effects are relevant when M > My /a, where 
My are the SM vector masses and a are their gauge couplings [10]. 

b) If (3<e<a ('lower triangle region') (3 is so low that its value is no longer relevant: R 
depends only on a/e showing, in the attractive case, a series of resonances. Indeed, by 
increasing a/e the potential develops more and more bound states: resonant enhancement 
happens when a bound state first appears with energy E B just below zero: in this case, 
at low energy R depends on K as dictated by a Breit-Wigner factor l/(K — Eb)- 

c) If e < < a ('upper triangle region') non-perturbative effects depend almost only on a/ (5: 
vector masses have a minor effect and R reduces to eq. (2). Its value does not depend on 
whether a loosely bound resonance is present or not. 

When DM DM bound states are relevant, the Sommerfeld correction might not encode the full 
dynamics: one should compute the density and life-time of DM DM bound states. 

We include one more effect, not discussed in [7]. At finite temperature T, the 7, Z, W 
masses and the mass splitting AM between the components of the DM multiplet are different 
than at zero temperature. First, because these masses are proportional to the SU(2)-breaking 
Higgs vev v, that depends on T and vanishes at T > T c when SU(2) invariance gets restored by 
thermal effects, via a second-order phase transition. This effect can be roughly approximated 
as [11] 

v(T) =vRe(l -T 2 /T 2 ) 1/2 . (3) 

The critical temperature T c depends on the unknown Higgs mass (T c ~ tuh for mu 3> v): we 
here assume T c = 200 GeV. 

Second, the squared masses of all SM vectors W, Z, 7, g get an extra contribution known as 
thermal mass. More precisely, even in the non abelian case, the Coulomb force gets screened 
by the thermal plasma: this can be described by a vector Debye mass, equal to [12] 

11 11 

m u(i) = Y 9y ' msu{2) = ~Q 92 ' msu(3) = 93 ^ 

in the SM at T ^> My/z- We approximate vector masses by summing the squared masses 
generated from v(T) with the SU(2)-invariant Debye squared masses of eq. (4). 

Let us now estimate the relevance of these thermal effects. For the W and Z bosons, thermal 
masses dominate at T > v, and for the photon and the gluons thermal masses dominate always. 
At temperature T DM particles have a typical energy K m T, so that thermal masses act like 
a contribution to e of order ej- ~ V^ra/? 2 - When this thermal effect is relevant, in fig. 1 one 
has to shift a/e vertically down to the diagonal of the shaded triangle: we see that (depending 
on the precise value of a/ (3) this shift has a mild or negligible effect on R, since thermal effects 
are relevant when R does not have a significant dependence on e. 

2.1 Computing Sommerfeld corrections 

Following [7, 10], we now give a operative summary about how to compute Sommerfeld cor- 
rections to s-wave DM DM annihilation rates. The set of two-body DMj DMj states that mix 
among them is labeled as i, j = {1, . . . N}. For example, when studying the wino triplet one 
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encounters the {DM~DM + , DM°DM } system. The dynamics is encoded in the N x N poten- 
tial matrix V and in another N x N matrix T, that describe the tree-level (co)-annihilation 
rates. The strategy is the same as in the abelian case of eq. (1), plus a careful book-keeping of 
indices. For each j one solves the following set of N coupled differential equations for tp^\r): 

i £)2„/,(i) N 

-j^- + EW = *W> (5) 

i'=l 

with boundary conditions: 

#(0) = <%, = y/M(K - V^cj~)4 j \oo) (6) 

chosen such that each wave has the same normalization at r = 0, where the annihilation occurs 
as described by the T matrices. The annihilation cross section <7j with given initial state i at r = 
oo is obtained by factorizing out the oscillating phase, = lim^oo ip\^ (r) / e l ^-\/ M< v K - v (°°)) r ^ 
and contracting with the annihilation matrix: 

a, = a(A ■ r • A% (7) 

The factor q is given by q = 2 if the initial state i has two equal DM particles (e.g. as in 
the DM°DM° state) and a = 1 otherwise (e.g. as in the DM DM + state). This formalism 
automatically takes into account that some states cannot exist as free particles at r ^> 1/Mw, 
when the kinetic energy K is below their mass. 

We next need to compute the Coulomb-like potential matrices V and the tree-level annihi- 
lation matrices V for all DM components. They split into sub-systems with given values of the 
L (orbital angular momentum), S (total spin) and Q (total electric charge) quantum numbers. 
Recalling that we neglect SM particle masses and that we consider s-wave (i.e. L — 0) tree-level 
annihilations into two SM particles, annihilations obey the following selection rules: 

0) two-body DM DM states with spin S = can annihilate into two SM vectors, that can 
have electric charge Q = 0, ±1, ±2. 

1) two-body DM DM states with spin S = 1 can annihilate into two SM fermions and into 
two SM Higgses, that can have electric charge Q = 0, ±1. 3 

The initial DM DM state with S = 1 of case 1) does not exist if DM is a scalar. 4 

We now give explicit expressions for the annihilation matrices and potential matrices 

Vuijji. Their off-diagonal entries do not have an intuitive physical meaning, and are precisely 



3 Non-perturbative corrections are mostly relevant when SU(2)l is unbroken: in this limit one could perform 
a simplified analysis by replacing the electric charge quantum number with the SU(2)l isospin quantum number 
/. Then, the only DM DM states that can annihilate into SM particles are those with S = and 1 = 1 and 
5 (that annihilate into SM vectors; 1=1 indicates the singlet) and the state with 5 = 1 and 1 = 3 (that 
annihilates into SM fermions and Higgses.). 

4 These selection rules also apply to DM DM annihilations relevant for indirect astrophysical DM signals, 
and reproduce well known results. E.g. if DM is a Majorana fermion the spin-statistics relation forbids S = 1, 
so that annihilations into vectors are not possible. A non vanishing amplitude for DM DM — > //, proportional 
to mf, is allowed when one takes into account the small mass rrif of the SM fermions /. 
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of, 



defined in terms of the imaginary part and of the real part of the two-body propagator ii 1 
where the indices ii' denote the two DM component in the initial state, and jf denote the two 
DM components in the final state. 

For MDM scalars with Y = one has (all states have S = 0) 



Vii';jf(r) 



33' 



8img D MM 2 A B 
{Mi + Mj - M)8 iV 5 jr + N ti ,N jf £ K ab (t£T*, + T A T B ) 
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where Nij — 1 if i ^ j and l/v^ if i = j. The indices A, B run over the SM vectors 
{7, Z,W + ,W~} with generators T A in the DM representation. The gauge couplings are in- 
cluded in the generators, such that the gauge-covariant derivative is = + iA A T A . We 
emphasize that the (off-diagonal) entries of T ± are defined up to an arbitrary phase normal- 
ization for each component of the multiplet: one needs to choose any convention such that the 
SU(2)-invariant DM mass term has the same sign for all DM components. The matrices Kab 
and Cab respectively describe SM vector propagation and the gauge content of all other SM 
particles, and are given by: 



Kab 



/l 0\ 
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1 

Vo 1 0/ 
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V 10/ 



41 2 
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For MDM fermions with Y = one has 



yS=0 



■pS=l 
ii',33 1 



Nii'Njji ST^rpArpB rl 

2^ 1 ii' 1 jj' L 'AB 



33 



87mg DM M 2 



A,B 



(9a) 
(9b) 



T s=1 is already summed over the 3 spin components. The potentials V remain the same as in the 
scalar case of eq. (8b), up to an extra — sign that must be taken into account when computing 
elements of V mediated by the W, such as Vo_,o- = — V0--0, because an extra — sign appears 
in the definition of two body fermion states: |DM°DM~) = — |DM~DM°). Furthermore, Fermi 
statistics forbids the existence of some states, such as |DM°DM°) with 5 = 1. 

Similar formulae apply to MDM candidates with 7^0, after taking into account that DM 
now becomes a complex particle (its normalization factor changes from iVoo = 1/2 to N q = 1) 
and that one is only interested in DM DM annihilations and two body states. Summing T^/ ^/ 
over all components one reproduces the total co-annihilation rates of [6], here reported and 
extended in appendix A. Explicit expressions for V and T are given in the next sections. 



3 Cosmological abundances 



To include non-perturbative corrections a dedicated lengthy analysis is needed for each MDM 
candidate. 
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Fermion triplet with ¥ = ('wino') M = 2.7 TeV 




DM mass in TeV Z = MIT 



Figure 2: The fermion triplet with zero hypercharge ('wino'). Fig. 2a (upper left): 
our result for its cosmological freeze-out DM abundance as function of the DM mass. Fig. 2c 
(upper right) show an example (for the indicated mass M) of the temperature dependence of the 
DM DM annihilation cross section, and fig. 2d (lower right) shows the resulting cosmological 
evolution of the DM abundance. Fig. 2b (lower left) shows our result for DM DM annihilation 
cross section relevant for indirect DM detection, as discussed in section 4- In each case, the 
continuous line is our full result, while the dashed line is the result obtained without including 
non-perturbative effects. 



The fermion triplet with Y = ('wino'). 

It is not automatically stable: one needs to impose a suitable symmetry. Since Y = the 
lightest component of this multiplet is neutral under the 7 and under the Z, and therefore it 
satisfies bounds from direct DM searches [13, 6]. This multiplet appears in supersymmetric 
models as SU(2)l gauginos, and behaves as MDM in limiting cases where it is much lighter 
than other sparticles. 

We include p-wave DM DM annihilations, but we only include non perturbative corrections 
to s-wave annihilations, that are splitted into four sectors with L = 0, total charge Q = {0, 1} 
and total spin S = {0, 1}. All sectors involve one component, except the sector with Q = and 
S = 0, that involves two components {DM~DM + , DM°DM }. In this case and in the following; 
we indicate the basis we employ by writing the charges of the DM components around the T 
and V matrices, and by indicating Q and S as pedices and apices on T and V. Using eq. (9), 
for the fermion triplet we get: 



ira 








- / 2A - A -\/2B\ 

Q=0 ~9M 2 \V2 2 y V Q=° ~ \-V2B )' 
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1 _25vra| yS=i = 2A _ A (n] 



Q=0 36M 2 

=0 na 2 n S=l 



V*:^ = A - B (12) 



rS^° 2 = S, ^r 2 ° = 2A + A (13) 



=1 9M 2' 0=1 36M 2' 

9 

0=2 - g^2- 

where A = 166 MeV, A = a cm /r + ct2C 2 N e~ Mzr /r and i? = a2&~ Mwr jr. These results are 
equivalent to those of [7]. 

We numerically solve the Schrodinger equations use the finite-temperature values for A and 
vector masses (e.g. A vanishes when SU(2)L-invariance is restored, etc.). Fig. 2a shows our 
result for its freeze-out cosmological abundance: as previously noticed in [7] non-perturbative 
corrections are relevant, and significantly increase the multi-TeV value of M that reproduces 
the measured cosmological abundance. In the supersymmetric case, a multi-TeV wino lightest 
super symmetric particle implies a fine-tuning in the Higgs mass above the 10 3 level, so that 
this scenario seems not motivated by the Higgs mass hierarchy problem. 

A bound state first appears for M > pa 2.5 TeV, in the DM DM system with total 
charge Q = and total spin 5 = 0. Fig. 2c and d show more details of the computation for 
M = 2.7 TeV. Since M is just above M*, the bound state is loosely bound, Eb ~ —67 MeV, 
and gives rise to a significant non-perturbative enhancement R of the annihilation cross section, 
apparent in fig. 2c as a dip at M = M*. The result is reliable, because temperatures T < \Eb\ 
(at which a dedicated treatment of bound states would be necessary) do not significantly affect 
the final cosmological abundance, as indicated by fig. 2d and by the fact that the dip is not the 
main effect. 



The scalar triplet with Y = 0. 

This MDM candidate is not automatically stable: one needs to impose a suitable symmetry. 
Having Y — 0, this candidate is compatible with bounds from direct DM searches. Any SU(2)l 
multiplet of scalars can have a quartic interaction with the Higgs, — \H(X^T a X)(W r a H) , that 
generates a tree- level mass splitting within the multiplet. Since this mass splitting is suppressed 
by the DM mass M [6] it is not unbelievable that scalars behave as MDM because Xh is small 
enough, \ H <C 0.05, that the A = 166 MeV mass splitting induced by gauge couplings is 
dominant. This is particularly plausible in the case of scalars with Y = 0, because A# is not 
generated by one-loop RGE corrections. Indeed RGE corrections that induce a A# proportional 
to g\ happen to vanish, because the Higgs if is a doublet, and the doublet representation of 
SU(2)l has vanishing symmetric tensor: Tr T^{Tj^,T^} = 0. 

To include non-perturbative corrections to DM DM annihilations we split them into three 
sectors, with L = 0, 5 = and total charge Q = {0, 1, 2}. The V and T matrices are 

+ +0 
5=0 _27m|-/3 V2\ S=0 _-(2A-A -V2B\ 

Q=0 ~9M 2 \y/2 2 y ^0=0-0 [-V2B )' li4J 

r g=o 2vra| V s=o = A + B r g=o 2vro| ~ = 

O-i gM 2 Q Q 9M 2 Q 

Fig. 3a shows our result for its freeze-out cosmological abundance: non-perturbative correc- 
tions are again important. Fig. 3c and d show details of the computation for M = 2.5 TeV, a 
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Scalar triplet with Y = M = 2.5 TeV 




DM mass in TeV Z = MIT 



Figure 3: The scalar triplet with zero hypercharge. Plots have the same meaning as in fig. 2. 



mass close to the critical mass M* « 2.5 TeV at which a bound state appears in the DM DM 
system with total charge Q = 0. 

We here do not study scalar and fermion triplets with |V| = 1; the neutral component has 
T3 = Y y£ and consequently couples to the Z, giving a too large direct detection rate. A 
non-minimal mixing with a singlet is needed to avoid this problem. One expects significant 
non-perturbative corrections, similarly to the Y = case, as SU(2)l gauge interactions are 
more significant than U(l)y the extra gauge interaction. 

The fermion quintuplet with Y — 0. 

This is the smallest multiplet whose lightest component is neutral under the 7 and under 
the Z and is automatically stable on the cosmological timescale. Indeed gauge couplings are 
the only renormalizable interactions (compatible with gauge and Lorentz invariance) that this 
multiplet can have with other SM particles. Its phenomenology is univocally dictated by a 
single parameter: its Majorana mass M. 

To include non-perturbative corrections to DM DM annihilations we split these scatterings 
into four L = sectors, with total spin S = {0, 1} and total charge Q = {0, 1, 2}. The Q = 
and S = system is composed by 3 states: 



++ /8A-4A -2B \ 

Vgr ° = + -2B 2A-A -3V2B , 
\ -3V2B / 



r S=0 = 3?ro 2 _ 

Q=0 25M 2 Q 
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Fennion quintuplet with Y = M = 10 TeV 




DM mass in TeV Z = MIT 



Figure 4: The fermion quintuplet with zero hypercharge. Plots have the same meaning as in 
fig. 2. 



Due to Fermi statistics, the Q — 0, S — 1 system has no DM°DM° state: 

++ + 

v s=i_ ++ (SA-4A -2B \ _tw4— f 4 2\ 

V Q=° ~ + I -2B 2A-A> i( 3=°~4M 2 - 2 IT [U) 



The systems with Q = 1 and 5 = {0, 1} have two states each: 



++ + 

t/ 5=o,i _ - ( 5A - 2A - V6B \ 
V Q =1 ~ I -y/EB A -35 J' 



++ + ++ + 

, 5=0 _ 3nal - / 6 v 7 ^ n s=i _ ^"i - ( 2 v 7 ^ 



=i 



25M 2 U6 1 ' Q=1 4M 2 U6 3 



Finally, in the systems with Q = 2 the tree-level annihilation rate T is non vanishing only for 
S = 0: 



v s =0 _0 ( 4A -2V3B] rS=0 _3vraio/ 4 ->/l2\ 



Q=2 +\-2V3B 2A + AJ' «= 2 25M 2 + ^-^12 3 

Fig. 4a shows our result for its freeze-out cosmological abundance: the value of M that re- 
produces the measured DM abundance was M ~ 4.4 TeV. It slightly decreases down to 
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Figure 5: Cosmological freeze-out abundance of the fermion doublet with Y = 1/2 ('Higgsino'). 
Plots have the same meaning as in fig. 2, except that since non-perturbative corrections negligibly 
affect the cosmological abundance we do not show details of the computations. 



M w 3.8 TeV after including p-wave and RGE corrections (perturbative result in fig. 4a), 
and increases up to almost M ~ 10 TeV after including the non-perturbative Sommerfeld cor- 
rections. Fig. 4d shows that Sommerfeld corrections are mostly relevant when T ^> AM, such 
that bound states are not expected to play a relevant role. 

The scalar quintuplet with Y — 0. 

As in the scalar triplet case, is plausible to assume that the quartic coupling with the Higgs 
is negligibly small, because it is not generated by RGE effects. The V and T matrices are 
related to those relevant for the fermion 5-plet as described in section 2.1. Non-perturbative 
corrections increase the value of the mass M that reproduces the cosmological abundance from 
M « 5.0 TeV [6] to M « 9.4 TeV. 

The scalar eptaplet with Y = 0. 

This is the smallest scalar multiplet that is automatically stable (its cubic scalar coupling 
identically vanishes) and contains a MDM candidate compatible with all bounds. Again, it is 
plausible to assume that the quartic coupling with the Higgs is negligibly small. Although we do 
not show dedicated plots nor the V and T matrices, we computed non-perturbative corrections 
finding that they increase the value of the mass M that reproduces the cosmological abundance 
from M ~ 8.5 TeV [6] to about M « 25 TeV. 

The fermion doublet with Y = 1/2 ('Higgsino'). 

It is not automatically stable, and, having Y ^ 0, is excluded by direct DM searches. Nev- 
ertheless, we consider it because n = 2 is the smallest multiplet, and because it is present in 
super symmetric models. The first problem can be solved by imposing a suitable symmetry, and 
incompatibility with direct detection experiments can be avoided by assuming a mass mixing 
with a neutral Majorana singlet (automatically present in supersymmetric models and known 
as 'bino'). We assume that this mass mixing is small enough not to affect the observables 
we compute. An almost-pure Higgsino LSP is realized in corners of the MSSM parameter 
space: it behaves as Minimal Dark Matter in the sense that its properties are not dictated by 
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super symmetry but only by gauge invariance. In this limit, the /x-term is the only relevant 
super symmetric parameter, and coincides with our DM mass parameter M. 

To include non-perturbative corrections to DM DM annihilations we split these scatterings 
into four L = sectors, with total charge Q = {0, 1} and total spin S = {0, 1}. The T and V 
matrices that describe each sector are 

+ 

s=0 Tiaj - ( 3 + 2t 2 + 1 4 3-2t 2 + t 4 \ 5=1 n<4 

1 n— n — _ . - ^ n .0 .A n . n . O , A i 1/ 



Q=0 64M 2 \ 3 - 2t 2 + t 4 3 + 2t 2 + t A J ' «=° 128M 2 

+ 




T/ s=o,i _ - | -a cm /r - (2c 2 - l) 2 e~ Mzr /Arc 2 -a 2 e~ M ^ /2r , 
K Q=° " 1 -a 2 e- Mwr /2r -e- Mzr a 2 /4rc 2 I 



s=0 _ Trnft 2 s=1 _ 25vra| S=0A _ a 2 e~ M ^ 2c 2 - 1 

iQ=1 "l6M 2 ' iQ=1 "64M 2 ' Kq=1 " + " r 4c 2 
where t = tan#w, c = cos^w and A = 341 MeV is the mass splitting among charged and 
neutral components generated by loop effects. Fig. 5a shows our result for its cosmological 
abundance: in agreement with previous studies (e.g. [6]) we see that the observed DM abun- 
dance is reproduced for M — 1 TeV. Non perturbative corrections are negligible. This also 
holds for the scalar doublet with \Y | = 1/2, so we do not study it. 

The Higgs potential in the Minimal Supersymmetric Standard Model depends on the \i 
parameter: a = 1 TeV gives a contribution to the squared Z-mass which is 2/i 2 /Mf = 240 
times too large: therefore the Higgsino as Minimal Dark Matter can only be realized at the 
price of a sizable fine tuning. 



4 Indirect detection of Minimal Dark Matter 

We now study the usual "indirect DM signals", generated by DM DM annihilations into 
SM particles. In the MDM case one only has tree- level annihilations into W + W~ , while the 
most interesting final states, 77 and 7Z (both detectable as 7 with energy practically equal 
to = M), arise at loop level. Assuming Y = the cross sections (averaged over DM 
polarizations) are, in the fermion DM case 

ff (DM° DM - W + W-)0 = (n 2 - l) 2 ^, a(DM° DM -> 77)^ = (n 2 - l?^^- 

W (23) 

and 2 times higher in the scalar DM case. These cross-sections can be significantly affected by 
non-perturbative Sommerfeld corrections [10], because in our galaxy DM has a non-relativistic 
velocity (3 ~ 10~ 3 . The computation employs the same tools already described in the compu- 
tation of DM DM annihilations relevant for cosmology. Non-perturbative effects can enhance 
the cross section by orders of magnitude for specific values of M: the values at which the 
potential develops a new bound state with Q,L, S = 0. In the triplet and quintuplet cases 
this enhancement is sizable, since the DM mass suggested by cosmology is close to one of such 
critical values. We show our results for the fermion triplet in fig. 2b, for the scalar triplet in 
fig. 3b, for the fermion quintuplet in fig. 4b, and for the fermion doublet in fig. 5b. The vertical 
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Figure 6: Spectrum 0/7 from DM DM annihilations around the galactic center, in a region of 
size Q = 10~ 3 , as seen by a detector with o~e/E = 0.15. The total rate is computed assuming 
the NFW profile, J = 1300, and for fermion M DM with n = 3 and M = 2.7TeV ; and for n = 5 
and M = 10 TeV. The continuous line is the H.E.S.S. result [14]. 



bands are the 3a range of M suggested by cosmology. We have not plotted annihilation cross 
sections into and ZZ since all MDM candidates with Y = predict 



These results allow to compute the energy spectra of annihilation products, such as photons, 
antiprotons, antideutrons and positrons. Note however that, precisely for the proximity to one 
of the critical values in mass, the cross sections depend strongly on the DM mass M and 
therefore the overall fluxes cannot be accurately predicted. Moreover, such total rates are 
affected by sizable astrophysical uncertainties. For example, the number of detected photons 
with energy E = M is [15] 



where is the angular size of the observed region and the quantity J ranges between a few and 
10 5 , depending on the unknown DM density profile in our galaxy. Keeping these limitations in 
mind, fig. 6 shows the predicted energy spectrum of detected photons from MDM annihilations. 
We have assumed realistic detector parameters (Q = 10 -3 and an energy resolution of 15%). The 
contribution from annihilations into W + W~ is included using the spectral functions computed 
in [15, 10]. The total rate has an uncertainty of about two orders of magnitude and is here 
fixed assuming the indicated value of M and the Navarro- Frenk- White DM density profile, 
J = 1300 [16, 1], which makes this signal at the level of the present sensitivity. Indeed the 
continuous line shows the spectrum of diffuse 7, dN^/dE ps 5 10 3 ( TeV/-E) 2,3 /m 2 yr TeV sr, as 
observed (in a region of angular size Q ~ 10~ 3 around the galactic center) by H.E.S.S [14] up 
to energies of about 10 TeV, where the number of events drops below one. Its spectral index 
suggests that the observed photons are generated by astrophysical mechanisms, rather than by 
DM DM annihilations. In view of the predicted value of cr 77 / o~ww '> the best MDM signal is the 
peak at = M. With the chosen parameters one gets iV 7 = 0.3 (2.6)/m 2 • yr for the fermion 
3-plet (5-plet). 



a^z = 2cr 77 / tan 2 #w = 6.5<x 



77' 



o~zz = tan 4 6\v = 10.8a, 



(24) 




(25) 
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5 Minimal Dark Matter at Ultra High Energies 



In this section we discuss the possibility that Minimal Dark Matter might be detected via 
the tracks left by its charged partners in experiments mainly devoted to Cosmic Ray (CR) 
detection, such as IceCUBE [18] and km3net [17]. Indeed MDM behaves differently from other 
neutral particles: 

• Neutrinos with energy E < 10 15 eV can cross the whole Earth [19] and, impinging on the 
rock or ice below the detector, produce a charged partner (a muon, a tau) that cruises 
the instrumented volume. 

• Neutral particles much heavier than neutrinos (e.g. the speculative neutralinos) interact 
rarely enough that they can cross the Earth even at Ultra High Energies (UHE), E 3> 
10 15 eV [20] . If they do interact before the end of their journey, they do produce a charged 
partner (e.g. the speculative chargino), that however usually decays almost immediately 
back to the neutral particle. For such reasons, it is generically difficult to detect them in 
this way. 

MDM is both heavy and quasi- degenerate with a charged partner, leading to the following 
behavior. To be quantitative, we focus on MDM multiplets with Y = 0. They have a stable 
neutral DM particle, DM , and a charged DM ± component which is 166 MeV heavier. The latter 
dominantly decays via DM ± -> DM°7r ± with a macroscopic life-time r = 44cm/(n 2 — 1) [6]. 
At Ultra High Energy Lorentz dilatation makes the DM ± life-time long enough that a sizable 
fraction of the travel is done as DM ± , leading to detectable charged tracks. 5 

This opens in principle an interesting and distinctive channel for the detection of DM. In 
order to assess its feasibility in practice, we need to study the system of MDM from production 
to detection. More precisely, first we study the issue of how a MDM system behaves when 
crossing the Earth at UHE, having assumed that the flux of UHE CR (observed up to E ~ 
10 20 GeV) does contain some DM particles. Mechanisms to produce such UHE DM are discussed 
later in section 5.2. We finally discuss detection signatures in section 5.3. 



5.1 Propagation of UHE Minimal Dark Matter 



Consider a flux of DM particles that is crossing the Earth at UHE. Via Charged Current 
(CC) interactions with nucleons of Earth's matter, DM ± particles are produced. Being charged 
particles traveling in a medium, these loose a part of their energy and eventually decay back 
to DM particles. This chain of production and decay is analogous to the process that tau 
neutrinos undergo in matter ("^-regeneration"). In the n = 5 case DM ±=t also play a role, but 
we will write explicit equations for the n = 3 (wino-like) case. Such a system is described by 
a pair of coupled integro-differential equations for the evolution with the position i of fluxes 
dNo/dE and n± = d(N + + N_)/dE of DM and (DM+ + DM"). They read: 6 



"o 



dn (e,E) 
di 



E 



T 



1 r°° 

n Q {E)a GC {E) + - / dE'n ± {E r 
I Je 



da cc (E',E) 



dE 



, (26a) 



5 Supersymmetric models with gravitino LSP may have a long-lived, charged next-to-lightest particle that 
behaves in a way qualitatively similar to MDM [21]. Furthermore, fermionic MDM with n = 3 and Y = arises 
in supersymmetry in the limit of pure-wino lightest supersymmetric particle. 

6 Similar equations have been of course used in the analogous problem for neutrinos [19] and for high energy 
neutralinos [20] . 
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dn±(e,E) M n±(E) M 

di = ~^^ + Nn 



\n ± (E)a cc (E) + jf" dE'n (E') E) 



+ 



where N^(i) is the number density profile of nucleons in the Earth [22], and E is the DM 
energy. The first terms account for DM 1 * 1 decays. The integro-differential terms within square 
brackets account for DM <-> DM ± CC scatterings. The last term accounts for the energy losses 
undergone by DM ± , due to electromagnetic and Neutral Current (NC) scatterings: 

dE\ oss _ dEi oss>em gLEi OSSj nc / 97 n 
di di di ' [ } 

We will later verify that the NC energy losses are subdominant and can be approximated as 
continuous, like electromagnetic energy losses. 

The cross section for the CC partonic scatterings DM d — > DM - u, DM°d — > DM + u, 
DM u -> DM + d and DM°w -> DM" d are all given by, for fermionic MDM: 



d&cc _ ( E da C c _\ 9l{n 2 - 1) 2M 4 + 2s 2 + 2si + g - 4M 2 g 
~~dT ~ \(s -M 2 ) d£ / 256vr (s-M 2 ) 2 ^-]^ 2 ,) 2 



(28) 



where s and t are the usual partonic Mandelstam variables. The total cross section for DM°iV — > 
DM+iV', DM-JV' is 

acc(s) = / dx q{x) / dt -j, (29) 

JO J-{s-M 2 ) 2 /s dt 

where q{x) is the parton distribution function [23], summed over all quarks and anti-quarks. We 
here approximated the neutron/proton fraction in Earth matter as unity. The nucleon N is at 
rest with mass m^, and DM has energy E , so that s = M 2 + 2m ^E and § = M 2 + 2xrri]yEo. 
A similar expression holds for NC scatterings of DM ± particles DM ± g — > DM ± g (whereas DM 
has no NC interactions due to Y = 0): 

da NC = 9t(gl q + 9 2 Rq ) 2M 4 + 2g 2 + 2st + g - 4M 2 s 
dt 16tt (s-M 2 ) 2 (t-M 2 ) 2 1 j 

where = T 3q — QgS 2 ^, Qr q = —QgS 2 ^ and q denotes any quark or anti-quark. The same cross 
sections, up to terms suppressed by powers of M\y,z/M, hold for scalar MDM. 

In fig. 7a we plot the MDM total CC cross section. We can understand its main features as 
follows. The total partonic cross section approximatively is 

9l(n 2 - 1) ( 2(xE/E*) 2 at xE < E* 
acc ~ 12MM 2 , ' \ 1 at xE » E* [6i) 

where E* = MM\y jm^ ~ 10 15 eV (the transition between the two regimes happens when 
t ~ — Myy)- Therefore one has a constant acc ~ 10 _33 cm 2 at xE ^> i?*. Upon integration 
over the parton densities, the nucleonic acc becomes a slowly growing function of E because 
at higher E partons with smaller x>E*/E contribute, and q(x) grows as x — > 0. This is 
analogous to what happens for neutrinos [24]. Indeed, at high energies, the MDM CC cross 
section in fig. 7 approaches the dotted red line of the analogous result for u^ T . 
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Figure 7: Left plot: CC cross section for fermionic MDM (black solid) with n = 3,5 on nucleons. 
The same cross section for neutrinos (green dot-dashed) is reported for comparison. Right plot: 
average energy loss parameters (3 = —dXnEjdl due to CC interactions (black solid lines), NC 
interactions (blue dashed), and electromagnetic interactions (red dotted). For comparison, the 
green dot-dashed line show the CC energy loss ofv^. The horizontal lines show the thickness 
of the Earth, crossed vertically. 



For the electromagnetic energy losses of DM ± one has the standard parameterization [25] 

dE\ oss em 



a + (3 em E 



(32) 



where a ~ 0.24TeV/kmwe approximates the energy loss due to ionization effects [25] and 
the j3 era term accounts for the radiative losses: it receives contributions from bremsstrahlung, 
e + e~ pair production and photonuclear scattering (i.e. inelastic electromagnetic collisions on 
nuclei). The latter two effects have comparable importance, while bremsstrahlung is subdomi- 
nant for a very heavy particle such as MDM [26]. We follow the approach of [26] and adopt a 
parameterization for (3 em with a mild dependance on the energy: 



ft 



2.5 10" 



TeV 



kmwe M 



1 + 



E 



10 15 eV 



0.215 



(33) 



Before proceeding to the numerical solution of the equations in (26), we can gain under- 
standing on the expected results by simplifying the treatment of the CC energy losses. In 
analogy with the standard parameterization, we can define an average energy loss suffered by 
DM particles due to CC interactions as 

- fcc(E) = = N N f dx a(x) f dt -i- % (34) 

E dl Jo J-(3-M*y/§ s-M 2 dt 

A completely analogous expression holds for /3nc — —1/E dE\ OSStNC /d£. Both are plotted in 
fig. 7b, in units of 1/kmwe = 10 _5 cm 2 /gram. The horizontal line indicates the value above 
which a particle looses a significant amount of energy while vertically crossing the Earth (its 
thickness being 1.1 10 5 kmwe). Energy losses can be understood analogously to the total cross 
section, and are dominated by partons with xE ~ E*. Parton distributions have been measured 
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at x > 10~ 4 and theoretical extrapolations to smaller x can have 0(1) uncertainties, that affect 
our results. 

The fraction of energy lost in one scattering is small, ~ < 10 -2 at the UHE energies 

that will be relevant for us, so that all energy losses can be approximated as continuous: 7 

- -jjj- = /Scc^o, — -jj- = a + (— + (3 cm + Pkc)E± (35) 

where Eq (E±) is the energy of DM (DM 1 * 1 ). Therefore, we do not need to study the evolution 
of the DM energy spectrum (dictated by eq. (26)) but we just need to follow how any initial 
energy decreases while crossing the Earth. This process is well approximated by the following 
system of two ordinary differential equations: 

— ? = -N N a cc N + 36a 

dl 1/r + l/N N a cc 

where < Nq < 1 is the fraction of DM present in neutral state, and 1 — N is the fraction of 
DM present in charged state. The initial conditions are Nq(0) = 1 and E(0) = E- m . 

Fig. 8a shows the numerical relevance of the various effects, and allows to understand 
their interplay. The DM interaction length starts to be smaller than the Earth thickness at 
E > 10 15 eV, but these interactions have no effects, since the produced DM 1 * 1 regenerates DM 
in a negligible length and with negligible energy losses. At E > 10 18 eV energy losses start to 
be moderately relevant, and at roughly the same energy the DM ^ mean free-path in Earth 
matter becomes comparable to the DM ± decay length: the Earth is crossed loosing an order 
one fraction of the initial energy, and spending an order one fraction of the path as DM ± rather 
than as DM . 

Fig. 8b shows the numerical results: the black continuous line shows the DM 1 * 1 fraction 
1 — iVo as function of the initial energy, after vertically crossing the Earth. The dashed line 
shows the same fraction as function of the final energy. Actually, we have two dashed lines 
obtained by solving the full and simplified system of evolution equations: their agreement 
confirms the validity of the continuum energy loss approximation. (In the full case we assumed 
an initial DM energy spectrum proportional to the CR energy spectrum). These DM 1 * 1 fraction 
is roughly given by 1 — iVo ~ (1 + Nncfcc/t)' 1 , and depends on the local nucleon density Nn, 
averaged over a typical scale of tens of km. This figure allows to compute the DM ± rate for 
any initial energy spectrum of DM. The upper red dot-dashed line shows the ratio E out /E- m 
between the final and initial energy. 

Fig. 9 shows the analogous result for the fermion 5-plet, which has a DM ±± component with 
a much shorter life-time ~ 0.05 cm. Due to its fast decay, the new DM ±± component is not 
present with a significant abundance at the UHE CR energies E < 10 20 eV. 

5.2 Production of UHE Minimal Dark Matter 

Minimal DM particles may exist at ultra high energies if they are a component of the primary 
cosmic rays, if they manage to be produced by standard UHE CR in the earth's atmosphere 



7 The fraction of energy lost in one DM — * DM decay is AM/M ~ 10 4 and can be neglected. 
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Figure 8: We consider the wino-like MDM: a fermion with n = 3 and M = 2.7 TeV. Left 
plot: The dotted lines show the length scales that characterize energy losses in matter with 
p = 6.3g/cm 3 . More precisely: with the upper blue dotted line we show the effective 1/Pcc 
experienced by a DM as if it never transformed in the charged state, and with the red dotted line 
the range of a stable DM ± . The solid lines show the effects of neutral to charge transformations 
and vice-versa: the solid ascending red line is the DM ± life-time while the solid descending pink 
line is the mean free path for DM — > DM ± CC interactions. Right plot: we consider a DM 
which crosses the Earth vertically with initial energy E in and we show E out /E in (upper red dot- 
dashed line) and the fraction of DM ± over the total number of DM particles fDM° + DM/* 1 ). 
These fractions are shown as a function of E ont in the simplified numerical approach (black 
dashed line) and in the full numerical approach (blue dashed line, almost superimposed to the 
former) as well as a function of E m (black solid line). 



or if more exotic CR scenarios will prove to be motivated. Let us look at these opportunities 
in turn. 

Although the mechanism that accelerates the UHE primary CR has not yet been estab- 
lished, the plausible astrophysical standard mechanism (first-order Fermi acceleration in mag- 
netic shock waves around various objects, such as gamma-ray bursts and supernova remnants) 
accelerates stable charged particles. In the MDM scenario, the neutral DM is accompanied 
by a charged partner with life-time r ~ cm: this value is not macroscopic enough to lead 
to acceleration of DM 1 * 1 particles. Moreover, the standard mechanism accelerates protons or 
charged nuclei in regions that are believed to be transparent: i.e. the local density is so small 
these primary CR particles negligibly hit on the surrounding DM accelerating it, or on the 
surrounding material producing DM pairs. 

Production of DM pairs might be relevant when UHE CR particles hit the Earth. The 
number of DM particles generated by one UHE proton with energy E> (2M) 2 /m p ~ 10 18 eV 
(which is the energy range that leads to a sizable DM ± fraction) is r ~ c/a pp ~ 10 -11 , where 
a pp ~ l/ m 7r i s the total hadronic cross section, and a is the pp — > DM DM cross section, 
computed in [6]. A UHE electron neutrino produces a higher fraction of DM, r ~ o ja v ~ 10 -6 , 
but only at higher energies E>(2M) 2 /m e . These processes lead to a double DM signature, 
but their rate looks too small for present detectors (see e.g. [20, 26]). 

A UHE DM rate detectable in forthcoming experiments can arise in more speculative sce- 
narios: e.g. they may be a component of the UHE CR generated by decays of ultra-heavy 
particle, in the so called 'top-down' scenarios [27]. 
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Figure 9: We consider the MDM fermion with n = 5 and M = lOTeV which crosses the earth 
vertically with initial energy E- m and final energy E out . We show E out / E- in (upper red dot-dashed 
line) and the fraction o/DM 1 * 1 and o/DM ±± over the total number of DM particles as function 
ofE out . 



In terms of absolute numbers, if DM constitutes a number fraction r of CR with energy 
above 10 17 eV (10 18 eV), one expects a flux of ~ 60r (« 6r) DM 1 * 1 per km 2 • yr, in the case 
n — 3. Atmospheric neutrinos generate a flux of ~ 10 3 (~ 10) up-going muons per km 2 - yr with 
E > 10 13 eV (10 14 eV): as discussed in the next section, this is a background to DM 1 * 1 searches. 

5.3 Detection of UHE DM ± 

Let us study how UHE DM 1 * 1 can be searched for. In a detector like IceCUBE or Antares DM ± 
roughly look like muons with fake energy = E±^cc/P^ ~ lO -4 .^, because muon energy 
losses are approximatively given by —dE^/dl = a + /3^E^ with (5^ « 0.2/kmwe [25]: indeed j3 
is roughly inversely proportional to the particle mass. 

One therefore needs to carefully study charged tracks to see the difference between a DM 1 * 1 
with E± ~ 10 18 eV and a muon with E^ ~ 10 14 eV. The muon would loose all its energy in about 
2 kmwe (with a characteristic energy loss profile) while DM" 11 have a much longer range in matter. 
On the contrary, the muon is essentially stable, while DM ± decays in about 10 km (less at lower 
energy: see fig. 8a), suddenly disappearing into a DM and a ^ with energy E n ~ m n E±/M. 
Furthermore, the characteristic scale for the DM «-> DM ± transformation is about tens of km. 
Unfortunately it seems very difficult to exploit these differences to discriminate between a ^ 
and a DM ± in the forthcoming IceCUBE km-size detector [18], because it only has a 1km size 
and modest granularity. 



6 Conclusions 

We considered DM that fills one SU(2) multiplet and only interacts with SM vectors ('Minimal 
Dark Matter' [6]). This model has one free parameter, the DM mass M, which is fixed from 
the observed DM abundance assuming that DM is a thermal relic. We have reconsidered this 
computation including non-perturbative corrections due to Coulomb-like forces mediated by SM 
vectors, that significantly increase M . For example, a fermion 5-plet automatically behaves as 
MDM (and in particular is automatically stable): its mass increases from M = 4.4 TeV [6] to 



20 




10 100 1000 10 4 

DM mass in GeV 



Figure 10: Predicted mass and predicted spin-independent cross sections per nucleon of MDM 
candidates. Scalar (F ermionic) SU(2) n-tuplets are denoted as n$ (np). The shaded region is 
excluded by [13] and the dashed lines indicate the planned sensitivity of future experiments [28]. 
The prediction for doublets with 7^0 hypercharge holds up to the caveats discussed in the text; 
furthermore we assumed the nuclear matrix element f = 1/3 and the Higgs mass m h = 115 GeV. 

10 TeV. The cosmological freeze-out abundance is plotted as function of M in fig. 4a (fermion 
5-plet), 2a (wino-like fermion 3-plet, M ~ 2.7 TeV [7]), 3a (scalar 3-plet, M ~ 2.5 TeV) and 5a 
(Higgsino-like fermion 2-plet, M « 1.0 TeV). We also discussed the scalar 5-plet (M pa 9.4 TeV) 
and 7-plet (M « 25 TeV). 

Having fixed M, we studied the MDM signals. 

Concerning the signals at colliders (see [6, 29]), we here point out that the charged compo- 
nents DM ± in the MDM multiplet manifest as non-relativistic cm-scale ionizing charged tracks 
negligibly bent by the B ~ Tesla magnetic field present around the interaction region. DM ± 
decays with 97.7% branching ratio into 7r ± , giving a relativistic track bent by the magnetic 
fields. Everything happens in the inner portion of the detector and this signature is free from 
SM backgrounds. However, at a hadron collider like LHC it seems not possible to trigger on this 
signature: level-1 triggers are located far from the interaction point and MDM candidates are 
so heavy that the event rate is too low. Unlike in other scenarios we cannot choose a favorable 
benchmark point in a vast parameter space. 

Astrophysics offers more promising detection prospects. 

The cross section for direct DM detection negligibly depends on M, and remains the same 
as in [6]: fig. 10 summarizes the situation. Since MDM is much heavier than a typical nucleus, 
direct DM searches can precisely reconstruct only the combination agi/M. 

Next, we considered indirect DM detection, computing how non-perturbative effects enhance 
the DM°DM° annihilations into W + W~, 77, 7Z, ZZ. Results are shown in fig. 2b (wino-like 
fermion 3-plet [10]), 3b (scalar 3-plet), 4b (5-plet), and 5b (Higgsino-like fermion 2-plet [10]). 
Signal rates can now be computed by combining this particle physics with (uncertain) astro- 
physics: for example fig. 6 shows the predicted spectrum of galactic 7. 

Finally, we assumed that some DM particles are present among the cosmic rays at Ultra 
High Energies and identified an unusual signal, characteristic of heavy and quasi-degenerate 
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multiplets containing neutral and charged components: a DM crosses the Earth without loosing 
much energy and, at E> 10 17 eV, a sizable fraction of the travel is done as DM 1 * 1 , leading to 
charged tracks, that might be detectable in neutrino telescopes such as IceCUBE or km3net [17]. 
Fig.s 8 and 9 quantitatively show how frequent this phenomenon is. 
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A Annihilation cross sections 



We here give results for the tree-level DM DM annihilation cross sections. We define the 
adimensional 'reduced cross section' 



/°*E^ (37) 
J-s ^ 8ns K ' 

where s, t are the Madelstam variables and the sum runs over all DM components and over all 
SM vectors, fermions and scalars, assuming that all SM masses are negligibly small. 
The DM abundance is computed by solving the Boltzmann equation 



where z = M/T, Ki is a Bessel function, Z = (1 — 1 , the entropy density of SM 

particles is s = 27r 2 ^ s T 3 /45, Y = n DM /s where hdm is the number density of DM particles 
plus anti-particles, and Y cq is the value that Y would have in thermal equilibrium. We can write 
a single equation for the total DM density because DM scatterings with SM particles maintain 
thermal equilibrium within and between the single components. We ignored the Bose-Einstein 
and Fermi-Dirac factors as they are negligible at the temperature T ~ M/26 relevant for DM 
freeze-out. 

We assume that DM fills one SU(2) multiplet with dimension n and hypercharge Y; when 
Y the conjugate multiplet is also added in order to allow for a gauge- invariant DM mass 
term. For example, the Higgsino has n = 2 and Y = 1/2; the wino has n = 3 and Y — 0. We 
define gx as the number of DM degrees of freedom: g% = n for scalar DM with Y = 0; gx = 2n 
for scalar DM with 7/0 and for spin- 1/2 (Majorana) DM with Y = 0; gx = for spin-1/2 
(Dirac) DM with F^0. For fermion DM we get 



9x 

° A 24vm 



(9C 2 - 21d)/3 + (lid - 5C 2 )P 3 - 3(2d(/3 2 - 2) + C 2 (/3 2 - l) 2 ) In i±| 



,' 3g 2 y - 1) + 20^F 2 g^-l)+Ag^ \( Q ^ 
~ 9x[ 16^ + 128^ ) V~T* (39) 



and for scalar DM we get 
9x 
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+9x 



20g*Y 2 



9$(n 2 - 1) 



48vr 



384tt 



/3 3 



(40) 



where x = s/M 2 and /? = wl — A/x is the DM velocity in the DM DM center-of-mass frame. 
The first line gives the contribution of annihilation into vectors, the second line contains the 
sum of the contributions of annihilations into SM fermions and vectors respectively. The gauge 
group factors are defined as 



c. 



J- TTT A T A T B T B = g ^ n yi + g 2£y 
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A.B 
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g* Y nY* + g\g 2 Y Y 



,n(n 



9\ 



9t 



n{n 2 — l) 2 
16 

n(n 2 — l)(n 2 — 5) 
16 



(41) 
(42) 



where the sum is over all SM vectors A = {Y, W 1 , W 2 , W 3 } with gauge coupling generators T A . 

The DM freeze-out abundance is accurately determined by the leading two terms of the 
expansion of for small j3, that describe the s-wave and the p-wave contributions. This approx- 
imation allows to analytically do the thermal average in eq. (38): 



<?A 



c s {3 + c p (3 3 + 



implies 



1A 



MT 3 e -2M/T 
32^3 



3T c s 
2M [Cp+ 2 



+ 



• (43) 



More complex models of DM can have extra interactions beyond gauge interactions, but the 
gauge contribution we computed accurately is the minimal model-independent contribution 
present in all cases. Eq. (38) also arises as a subset of the Boltzmann equations that describe 
leptogenesis from decays of gauge-interacting particles: the scalar triplet corresponds to n — 3, 
Y = 1 [30], and the fermion triplet to n = 3 and Y = [31]. 
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